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ABSTRACT 

Transverse stratification is a common intrinsic feature of astrophysical jets. There is growing evi- 
dence that jets in radio galaxies consist of a fast low density outflow at the jet axis, surrounded by a 
slower, denser, extended jet. The inner and outer jet components then have a different origin and launch- 
ing mechanism, making their effective inertia, magnetization, associated energy flux and angular momen- 
tum content different as well. Their interface will develop differential rotation, where disruptions may 
occur We here investigate the stability of rotating, two-component relativistic outflows typical for jets 
in radio galaxies. For this purpose, we parametrically explore the long term evolution of a transverse 
cross-section of radially stratified jets numerically, extending our previous study where a single, purely 
hydrodynamic evolution was considered. We include cases with poloidally magnetized jet components, 
covering hydro and magnetohydrodynamic models. With grid-adaptive relativistic magnetohydrodynamic 
simulations, augmented with approximate linear stability analysis, we revisit the interaction between the 
two jet components. We study the influence of dynamically important poloidal magnetic fields, with vary- 
ing contributions of the inner component jet to the total kinetic energy flux of the jet, on their non-linear 
azimuthal stability. We demonstrate that two-component jets with high kinetic energy flux, and an inner 
jet effective inertia which is higher than the outer jet effective inertia are subject to the development of a 
relativistically enhanced, rotation-induced Rayleigh-Taylor type instability. This instability plays a major 
role in decelerating the inner jet and the overall jet decollimation. This novel deceleration scenario can 
partly explain the radio source dichotomy, relating it directly to the efficiency of the central engine in 
launching the inner jet component. The FRII/FRI transition could then occur when the relative kinetic 
energy flux of the inner to the outer jet grows beyond a certain treshold. 

Subject headings: ISM: jets and outflows - Galaxies: jets - methods: numerical, relativity 



1. Introduction 

There is strong observational and theoretical evi- 
dence that magnetic fields play a crucial role in the 
acceleration and the coUimation of extragalactic jets. 
Most Active Galactic Nuclei (AGN) jet formation 
scenarios involve magnetic fields threading a rotat- 
ing black hole (in the ergosphere) and its accretion 
disk, thereby removing from them angular momen- 
tum, allowing the central black hole to accrete. At 
least close to the jet launching region, the jet rota- 
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tion profile persists, reflecting its (general relativistic 
and/or magneto-rotational) origin. Thus both ingredi- 
ents, magnetic fields and rotation, are very important 
in jet formation, as well as in jet propagation and sta- 
bihty 

Moreover, detailed astrophysical jet observations 
point out that relativistic jets are structured, in a di- 
rection perpendicular to the jet axis, typically con- 
sisting of a fast spine and slower outer flow. In the 
case of AGNs, this jet structuiing plays an important 
role in explaining the morphology of the jet high en- 
ergy radiation ("GhiselUni et al.'2005',!Hardcastle|2006 



Jester et al. 2006 2007; Siemiginowska et al. 



2007 



Kataoka et al.||2008| l, with sometimes clear evidence 



for a very fast, light inner jet and a heavy slow outer 
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outflow ( |Giroletti et al.||2004[ l. Furthermore, obser- 
vations of the TeV BL Lacertae objects show bright- 
enings and rapid variabiHty in their TeV emission. 
This variation in their high energy emission impHes 
high Lorentz factor flows occuring at smafler scale, 
suggesting uhra-relativistic bulk motion of the (inner) 
jet. At the same time, complementary (radio) obser- 
vations with VLBI of the pc-scale jet structure indi- 
cate a broad, 'slowly' (albeit relativistic) moving out- 
flow. In combination, this clearly suggests the pres- 



ence of a two component jet morphology (Ghisellini 



et al.|[2005] ). Two-component jet structure has also 
been proposed in more theoretical work, addressing 
the physics of jet launching, collimation and propaga- 
tion mechanisms (Bogovalov & TsinganospOOl Sol 



|etal.|1989t|Meier|2003| l 

While our jet dynamics computations will be rep- 
resentative for AGN jet conditions, radially struc- 
tured jet flows are now known to exist in virtually 
all astrophysical jet outflows. Transversely structured, 
ultra-relativistic jet-like outflow has been proposed 



in the context of Gamma Ray Bursts (Racusin et al. 
[2008'), to explain the break observed in their after- 
glow light curve. In the case of stellar outflows, re- 



cent observations of some T Tauri jets (Bacciotti et 



aLl [2000t [Gtinther et HI][2009 ^ also suggest a fast in- 



ner outflow bounded by a slow outer outflow. In these 
young stellar objects, a clear signature of jet rotation 
around the symmetry axis was detected (Bacciotti etl 



al.|2002[[Wortas et al.|2005[ [Coffey et al. 2004, 2007; 



fully supporting scenarios of magneto-centrifugal jet 
launch and acceleration. Theoretical models of two- 
component jets in classical T Tauri ([Bogovalov & 



Tsinganos 2001] [Mehani et al.|2006b[ |Cranmer|2008 



Fendt|2009| ) then postulate that the inner outflow is tur- 
bulent and pressure driven, associated with the young 
star wind. The inner jet then has a small opening an- 
gle, as it is collimated by the outer jet, which is in 
turn magneto-centrifugally driven from the surround- 
ing disk. The outer disk wind then carries most of the 
mass loss in the jet. Various authors (Meliani et al. 



2006b Fendt|2009 1 have demonstrated using axisym- 



metric MHD simulations that the outer outflow is self 
collimated by its intrinsic magnetic field, and that the 
turbulent inner outflow gets collimated by the outer jet. 



Furthermore, Matsakos et al. ( 2008 1 investigated the 
topological stability of two-component outflows for 
young stellar objects, performing extensive numerical 
simulations to determine whether analytic self-similar 
models demonstrate robustness in axisymmetric con- 



ditions. 

Also for relativistic jet simulations, axisymmetry 
assumptions are often adopted, excluding the devel- 
opment of ah non-axisymmetric perturbations. These 
can address details of how helical field configurations 
(naturally expected from magneto-centrifugal launch 
mechanisms) effectively may transport their helicity 
down the jet beam (Keppens et al. 2008), with magnet- 



ically aided reacceleration by field compression across 
internal cross-shocks. While Keppens et al. (20081 



concentrated on kinetic energy dominated jets, ini- 
tially Poynting flux dominated jets were simulated by 



Komissarov et al. ( 2007 1 in axisymmetric relativistic 



MHD, finding that the transition to a matter-dominated 
jet regime occurs very close to the central engine 
(within O.Olpc). Our model computations will there- 
fore assume kinetic energy dominated jets. 

As far as the magnetic field topology is concerned, 
we will restrict ourselves in this paper to purely 
poloidally magnetized jet components. Our two- 
component jet model determining our initial condi- 
tions can actually allow for helical fields, as explained 
in Sect. [2] (we include this more general case here, for 
future reference in follow-up studies). As indicated 
before, during the first acceleration phase of AGN jets, 
magneto-centrifugal mechanisms play an important 
role and a helical or even strongly toroidal magnetic 



field is likely produced ( Fendt |1997[ [Meliani et al. 
12006a i McKinney & Blandford 2009; Komissarov et 



ial..2007) . ,McKinney & Blandford, ( ,2009) present 3D 
general relativistic MHD simulations for rapidly ro- 
tating black holes producing jets with strong toroidal 
fields. They find a prominent role of the accreted mag- 
netic field geometry for achieving 'stable' jets. Helical 
or strongly toroidal field topologies can be subject 



to current-driven kink instabilities (Begelman 1998 1, 
with m = 1 toroidal modes that helically displace 
the jet axis. This requires full 3D numerical simula- 
tions, such as performed by |Baty & Keppens ( 2002|l in 
non-relativistic MHD, or addressed by [Mizuno et al. 



( [2009 ^ in relativistic MHD for a static force-free equi- 
hbrium. Dispersion relations for non-axisymmetric 
modes and m - \ kinks in particular for relativistic 
MHD were analysed by 'Begelman] ( [1998| l for purely 
toroidal fields, and electromagnetically dominated. 



force-free jets were analysed spectrally by Istomin 



& Pariev ( 1996) 1 and more recently by Narayan et al. 
P009» . 



In our work, we will restrict attention to 2.5D sce- 
narios with the somewhat unusual assumption of trans- 
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Fig. 1. — 3D schematic view of the overall AGN disk- 
jet configuration (indicating the accretion disk and the 
two-component jet). We model the jet evolution in the 
transverse plane. 



lational invariance along the jet axis. The overall con- 
figuration is schematically indicated in Fig. [l] and we 
simulate a transverse cross-section of the jet at suf- 
ficient distance from the two-component jet source, 
where all three velocity components (axial, azimuthal 
and radial) are included, but their variation along the 
jet axis is ignored. Our aim is to investigate all non- 
axisymmetric instabilities, primarily induced by the 
(sheared) rotation. This approximation is valid be- 
cause in the poloidal direction, the flow is supersonic 
with high Lorentz factor, and then the growth rate of 
poloidal instabilities is expected to be low. On the 
other hand, the rotation is subsonic, facilitating the 
growth of toroidal instabilities. We then address jet 
stability in a cross section of a rotating two-component 
jet, initially collimated by thermal pressure and/or 
poloidal magnetic field. Consecutive snapshots of the 
cross-sectional evolution can be interpreted as mim- 
icking the jet flow conditions at increasing distance 
from the source. Note that this particular assumption 
allows us to follow both axisymmetric and all non- 
axisymmetric (including m - 1) mode development, 
but does exclude helical mode axis displacement typ- 
ical for kink modes. Our model therefore mimics jet 
evolution adequately as long as the radial axis dis- 
placement is smaller than the axial wavelength associ- 
ated with possible m - \ kinks. Our assumption does 
also neglect conical jet expansion, assuming cylindri- 
cal propagation. This is justified given the low ob- 



served values for jet opening angles. 

Since we defer the study of toroidal and helical 
magnetic field configurations in 2.5D and 3D to later 
work, we start off with numerically investigating the 
influence of a purely poloidal magnetic field on the sta- 
bility of rotating, two-component relativistic jets. Our 
work complements the studies looking into kink de- 
velopment, by putting the emphasis on the azimuthal 
variation and on the eff'ect of the two-component jet 
stratification. As far as a purely poloidal magnetic field 
topology is concerned, the work by'Spruit et al. ( 1997 1 
suggests that jets should be collimated by poloidal 
magnetic field pressure, rather than by toroidal mag- 
netic field, as toroidal jet magnetic fields can intro- 
duce kink instability (but only slow mode growth was 



found for force-free jets by Narayan et al. ( 2009 1). To 
justify purely toroidal fields, we can argue that the 
toroidal magnetic field in the jet acceleration phase 
can be gradually dissipated involving reconnection, a 
mechanism which in turn contributes to (axial) jet ac- 
celeration ( Spruit 2008, |Meliani et al.|2006bl l. Indeed, 
during the acceleration phase a fraction of the Poynt- 
ing flux (angular momentum) carried by the magnetic 
field is converted to kinetic energy by internal dissipa- 



tion/reconnection of the toroidal magnetic field ( Spruit 
|2008[|Mehani et al.|2006b| . In accord with this mech 
anism, we model the region where the eventual jet ro- 
tation is low. 

In this paper, we analyse five cases in detail, to de- 
termine the effects of differing poloidal magnetic field 
configurations in the two-component structure on its 
long-term stability, and on the overall deceleration of 
the jet as it propagates away from the central source 
regions. The role of poloidal magnetic field is in these 
cases most prominent in its added effect on total pres- 
sure and effective fluid inertia, and this is shown to 
play a prominent role in the two-component jet stabil- 
ity. 

2. Two-component jet model 



We set up a relativistic magnetohydrodynamic 
model of a two-component jet, which elaborates on the 
earlier model described in Meliani & Keppens ( 2007| l. 
The jet model uses as its basic characterising parame- 
ter a total jet kinetic luminosity flux typical for pow- 
erful radio loud galaxies, namely Ljet,Kin = lO'^^ergs/s 
( [RawUngs & Saunders |[T99T] |Tavecchio efaLjllOO? i. 
Secondly, the outer radius of the two-component jet 
is taken to be Rq^i ~ O.lpc. This is a value directly 
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inferred from observed values for M87 ( [Biretta et al.| 
|2002 |), which has a known jet opening angle of 6° at 
Ipc from the jet source. For the less constrained in- 
ner jet radius, we adopt the initial value - ^out/3. 
The remainder of the t - Q condition is then character- 
ized by the velocity profiles in both inner and outer jet 
components and by the magnetic configuration in the 
two-component jet. 

2.1. Initial flow profiles and Rayleigh criterion 

For the initial flow condition, we adopt a uniform 
outflow velocity along the jet axis in each jet com- 
ponent, with Vz^out - 0.9428 (hence yz ^ut ~ 3, since 
speeds are normalised to the light speed) for the outer, 
slow jet located within Ri^ < R < Rom- This represents 
a typical outflow value for relativistic jets at parsec 
scale from the source. The inner jet has a much faster 
outflow speed, which we set to Vz,in = 0.99939, with 
corresponding Lorentz factor yz.in ~ 30. The differ- 
ence in Lorentz factor by an order of magnitude in the 
inner versus outer jet layers is representative (though 
somewhat exaggerated) for differences inferred for the 
BL Lac object Markarian 501 between its central spine 
and surrounding (shear) layer jet, as described by |Giro-| 
|letti et alT] ( |2004[ ). However, it remains lower than the 
Lorentz factor y ~ 50 suggested for TeV blazar PKS 
2155-304 dGiroletti & Tavecchio|2008) . 

Two-component jet models also commonly assume 
that the spin of the inner beam is higher than the spin 
of the outer jet. The physical argument is that the in- 
ner beam is extracting angular momentum from the 
inner parts of the accretion disk and from the black 
hole itself, both implying fast rotation. Therefore, we 
adopt different initial rotation profiles for inner ver- 
sus outer jet, both radially self-similar. Although we 
will in our numerical results below only adopt purely 
poloidal magnetic field configurations, i.e. we will 
have B = B^{R, <p; t)i^, one can generalize our initial 
equilibrium configuration to allow for both toroidal 
speeds and toroidal magnetic fields, sharing the same 
/^-dependence. In that general case, one takes in prac- 
tice 



— 



b ■ f^V 



am/2 
fl„,„/2 



R<Rin, 
Rin < R < R„,,t , 



R<Rh 

Ri„ < R < Rm. 



(1) 



(2) 



where the V = V^e^ + VzS, is the 3-velocity, and sim- 
ilarly for the magnetic field. At f = 0, we then have 
a discontinuity in the toroidal velocity at the bound- 
ary between the two jets when fixing v^_i„ = 0.01 
and v^out - 0.001. A faster rotating inner compo- 
nent is consistent with a small expansion of the inner 
streamlines, which are confined by the outer outflow. 
The fast rotating inner jet effectively extracts the an- 
gular momentum from the central region and carries it 
away with a very small mass flux. In fact, our mod- 
els typically have the inner jet extracting order 10% 
of the total angular momentum associated with the 
two-component jet. We fix the exponent ai„ - 0.5 
for the inner jet and = -2 for the outer jet. We 
then choose for the outer jet a rotation profile im- 
plying a radially constant angular momentum. This 
is motivated by the fact that this component is be- 
lieved to be launched from an accretion disk, where 
the outer streamlines in the jet expand faster than the 
inner streamlines. 

The rotation profile of the inner jet then obeys the 
relativistic equivalent of the Rayleigh criterion for sta- 
bility, namely that the angular momentum flux given 
by 



7- -^V^R- 



Bp 

R B,n . 



(3) 



P ry "p 

increases with R ( Meliani et al.||2006a[ l. This expres- 
sion uses proper density and pressure p, p, assumes 
a polytropic index F appearing in a simple polytropic 
equation of state (our actual numerical simulations will 
in fact relax this assumption and use a full Synge-type 
equation of state), and writes Vp, Bp for poloidal three- 
velocity and magnetic field strengths. For the chosen 
radial prescriptions for given by Eq. ([TJi, we have 
for the inner jet 



dR 



(fl,„/2+l)>0, 



(4) 



making the inner jet centrifugally stable. The outer 
outflow on the other hand is marginally stable, follow- 
ing the same argument. In the numerical simulations 
further in this paper, we neglect toroidal field B^, and 
then the interface between the two components does 
not verify the Rayleigh criterion and this shear flow 
interface is unstable from the start, since the angular 
momentum at the interface decreases when 



(Pout + p^^"'" 1 /'out) (pin + p.J^i/'in) 

Tout i-V,out < Tin ' 



Pout 



(5) 
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This means that in all cases we will simulate, the initial 
shear flow interface at = Ri^ will develop small-scale 
instabilities, and we intend to address the ultimate non- 
linear stability of such initial two-component structure. 
The central question investigated then is whether the 
dynamics will lead to complete destruction of the ini- 
tial two-component nature of the outflow at further dis- 
tances from the source. Related to this question, we 
will quantify the potential deceleration of the central 
fast jet by mixing processes, induced by the nonlinear 
evolution of the two-component jet. 

2.2. Initial densities and magnetic configuration 

The densities and the poloidal magnetic field are 
assumed to be constant within the three different re- 
gions, namely throughout inner jet beam, outer jet, and 
external region. The external medium is hot and rari- 
fied and has a small reference number density. The 
typical ISM number density value of Icm""^ is in the 
computation used as a scaling value. We arbitrarily set 
Pmed = 'WpWrned with «n,ed - lO^^cm""', and adopt c = 1 
and a unit of length of Ipc. This sets the unit of pres- 
sure, and a corresponding scaling value for the mag- 
netic field is then 0.1375 Gauss (as mass is normalized 
to proton mass nip). 

The jet environment represents a relativistically hot, 
rarified, static external medium. This is consistent with 
the fact that at the parsec scale, the jet head has pre- 
viously shocked the jet surroundings, so that the jet 
itself gets embedded in a dilute, hot medium. As- 
suming a number density for this external medium of 
«med = lO^^cm"^ is appropriate, since this is near val- 
ues obtained in numerical simulations of jet propaga- 
tion and are then typical for the jet cocoon surround- 
ings ( Mehani et al.||2008 l. Values for the jet compo- 
nent proper densities are estimated from the kinetic lu- 
minosity flux and its division over inner and outer jet 
components. Requiring that the outer jet carries a case- 
specific fraction of the total kinetic luminosity flux by 



'jet,Kin 



(Tout /lout - l)Poutrout7r(^out - ^?n)Ki 



(6) 

we can deduce the density pout, if we in addition pre- 
scribe the (outer) jet Lorentz factor and make the initial 
approximation that the thermal energy in the outer jet 
is negligible compared to its mass energy (as valid in 
all cases studied). In Eq. (j6]), /lout is the specific en- 
thalpy in the outer jet, and its expression depends on 
the equation of state. When we consider that the in- 
ner component carries only a small fraction = 



of the total kinetic luminosity flux, such that the outer 
jet carries the remaining /out = 99%, we find that the 
resulting density is pout ~ 119.94 x lO^med- A sim- 
ilar argument for the inner, light jet, works out to fix 

Pin ~ 6.92p,ned- 

While in all cases we investigate here numerically, 
we assume these initial, piecewise constant, densities 
for the inner and outer jet, the actual kinetic luminos- 
ity carried by each component will change from one 
case to another. This relates to our prescription of the 
pressure and magnetic variation in the jet components, 
and we will simulate five cases with varying thermal 
pressure contribution at the jet axis. In four cases (A), 
(B2), (C) and (D), the total pressure at the jet axis ends 
up similar, while case (Bl) has a lower total pressure 
value. By varying the relative contribution of thermal 
to magnetic pressure in the jet components, we aim 
to analyze the effect on overall jet stability of varying 
magnetic field and kinetic energy contributions to the 
kinetic luminosity in each component. In all five cases 
studied, we set the pressure of the external medium to 
ensure total pressure balance with the outer jet. Also, 
in the four cases (A), (B 1 ), (C), (D), we assume that the 
thermal energy in the inner (hot) component is higher 
than its mass energy. Only in the case (B2), the ther- 
mal energy is lower than the mass energy. They differ 
in the following sense: in the first case (A), both com- 
ponents are non-magnetized and the pressure at the jet 
axis is po = 2.3 in our simulation units. This purely 
hydrodynamic case differs from the hydro case already 
studied in Meliani & Keppens ( 2007| , as now the in- 
ner jet has an even higher Lorentz factor, a slightly 
different inner rotation profile, and an on axis pres- 
sure value which makes the entire inner jet relativis- 
tically hot. In the cases (Bl) and (B2), only the in- 
ner jet is magnetized and we set it to a constant value 

Bz.in = ^OXHy^/^i^ ~ 0.789 (making the jet param- 
eter cr - Blliy^p) expressing magnetic to rest-mass 
energy of order 0.01, as for kinetic energy dominated 
jets) in the case (Bl) and B^jn ~ 2.28 in the case (B2). 
In the case (Bl), the thermal pressure on axis is set to 
po = 0.1 and in the case (B2) the thermal pressure on 
axis is po - lO""*. In the case (C), only the outer jet 

is magnetized with Bz.out = -^0.005yoy(Pout ~ 2.323, 
and we adopt the same order on axis inner pressure 
as in (A), namely po « 2.3. In case (D), both inner 
and outer jet are magnetized, with in ~ 0.789 and 
fiz.out ~ 2.323. We then take the pressure at the jet 
axis pd ^ 2. In all magnetized cases, the magnetic 
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field strenght is set according to the observations at 
parsec-scale ( |0' Sullivan & Gabuzda|[2009l ). Clearly, 
in all four cases (A), (Bl), (C) and (D), the thermal 
energy in the inner component dominates, only in the 
case (B2) the thermal energy is lower than the mass en- 
ergy. The pressure profile through inner and outer jet 
is taken from transverse equilibrium conditions, dis- 
cussed next. 

2.3. Two-component jet pressure profiles 

Since we assume an initial near steady-state ax- 
isymmetric two component jet, the actual pressure 
variation is deduced from transverse equilibrium be- 
tween pressure gradient, centrifugal force, and Lorentz 
force. This is expressed by the steady-state equation 



( [MeUani et al.|2006a[ [^pl & Camenzind|1993[|Hey^ 
Ivaerts & Norman|2003| l 



(7) 



Here p is the proper density of particles, p and h are the 
pressure and enthalpy per particle, respectively, is 
the fluid three-velocity, y the Lorentz factor, and (E, B) 
denote the electromagnetic fields. We have pe = V-ii 
and for steady-state conditions = V A B as the asso- 
ciated charge and current densities, while E - -V A B 
in ideal relativistic MHD. 

If we introduce the total pressure in the fluid frame, 



B^-E^ 



= P + 



Bl + Bl {B,V,-V,B,f 



the radial component of Eq. (|7]l writes as 



(8) 



di// 
dR 



r y% 



2t/2 



r- 1 " R 

.2t/2 \ 



r Bl 



1 + 



2(r-i)j 



r- 1 2 
{-bI + {b,v,-v,b.)^). 

(9) 



This has assumed a polytropic equation of state. We 
can then determine the radial pressure profile p{R) 
by solving the first order differential equation (j9]l and 
using ([8]l, together with the self-similar prescriptions 
from Eq. ([TJ-Q to get 



Ri, 



(r-1) 



(r- l)(a-H2) / 2 /, 



Ru 



(10) 



In this expression, the constant a is a self-similarity 
exponent, which we use to set the inner jet exponents 



flin and the outer aout, and the constant a ■ 



i-V: 



, which 



is different for inner versus outer jet a\n, ffout- Also, f 
is an integration constant, which is to be deduced from 
boundary conditions. For the inner jet, we use it to set 
the pressure on axis po from 



(Fin - 1) 

^in = PQ+ Pin 

A in 

. (r-i)(fl + 2) 



2a{a (r- i) + r) 



(11) 



To fix the constants for the outer jet, we have to match 
conditions at the interface between the two compo- 
nents. At the contact between the two components, 
Rankine-Hugoniot conditions should hold. In terms of 
effective total pressure ([8]l obeying Eq. (|9]l, we need to 
match at the interface 



tAin(^in) = lAout(^in) ■ 



(12) 



This yields for the constant pressure integration in 



Eq. (10 1 for the outer jet 

Tout ~ 1 



l/'in(^in) + 



"Pout ~ 



out ^ 



"^^,OUt (^(^ ,out ^z.out ~ ^(^.OUt^Z.OUt) 



(l-aout)'-(^»'-') 

21 



^ '!_ (1 - gput) (rout - 1) (flout + 2) 
\ 2 2 QTout [flout (Font ~ 1) + Tout] 



(1 -C^out) 

(13) 



This fixes the pressure variation throughout both com- 
ponents. In the external medium, the thermal pressure 
is constant and set equal to the total pressure in the 
outer jet component. This is actually a slightly dif- 
ferent value for cases (A), (Bl) and (C), (D); albeit 
of the same order. The governing equation of state is 



taken as a Synge-type relation, also used in Meliani et 



al. (2004 1. As a result of the above profile prescrip- 



tions, the outer jet is relativistically cold with an ef- 
fective polytropic index Fgff - 5/3 and the inner jet 
has ultra-relativistic state Feff = 4/3. Case (B2) is an 
MHD variant to the hydro case with more stratified 
inner effective polytropic index already simulated in 
|Mefiani & Keppens] ( |2007[ ). The matter state between 
the two components is different (relativistically hot in- 
ner and cold outer jet) and we will see that during the 
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time evolution, a shear region with intermediate mat- 
ter state forms. This variation in the matter state be- 
tween various jet regions makes it vital to use a Synge- 
EOS to model growing instabilities and resulting tur- 
bulence. However, in the initial conditions of all cases, 
the polytropic index is relatively constant throughout 
each component. This is why we could use the poly- 
tropic EOS assumption to deduce a near-equilibrium 
solution of the initial two-component jet. 

2.4. Dimensionless characterization and numeri- 
cal setup 

As a result of the initializations described above, 
the local fast magnetosonic speed in the inner jet for 
case (D) and the local sound speed in the inner jet for 
cases (A) and (C) are of the order of 0.6 (light speed 
units), while cases (Bl) and (B2) have a higher local 
fast magnetosonic speed. In the outer jet, the local fast 
magnetosonic speed of cases (C), (D) and the local 
sound speed in case (A) and (B2) are of the order of 
0.2, while it is 0.07 in the case (Bl). Both inner and 
outer jets are kinetically dominated. The inner jet has 

y V 

an effective relativistic Mach number Mfast - 
up to ~ 40 for cases (A), (C), (D), while the outer 
jet has Mfast ~ 14.0 in cases (A), (B2), (C), (D) and 
■Wfast ~ 35 in case (Bl). However, both components 
are subsonically rotating, and during a complete rota- 
tion of the inner jet, the fast magnetosonic or sound 
wave will propagate about 200 times from the edge of 
the inner jet to the axis and back. Such overall config- 
uration can easily develop non-axisymmetric instabil- 
ities, with growth times of order of the radial sound- 
crossing time ( |Hardee|2004) l. 

The computational domain of this simulation is a 
2D box of size -0.3pc < x < 0.3pc and -0.3pc < 
y < 0.3pc. The simulation is performed in Cartesian 



coordinates using an HLLC flux formula (Mignone & 
|Bodo|[2005] l. HLLC is a two state extension of the 
Harten, Lax and van Leer flux formulation (HLL), 
which includes a proper representation for the con- 
tact wave. We use a piecewise parabolic method 



(PPM) limiter ( [Mignone et al.||2005[ ). The combina- 
tion of PPM reconstruction (third order accurate) and 
HLLC flux computations is extremely robust and han- 
dles both sharp discontinuities and turbulence develop- 
ment accurately. The lateral boundaries assume open 
boundary conditions, with a clipping of any inwardly 
directed momentum as soon as turbulent flow features 
start crossing the boundaries. The simulation is run 



till time t - 50, which due to our normalization trans- 
lates to 163.2 year The corresponding distance of jet 
propagation of the jet beam during this time is about 
50pc. The simulation is done using the AMRVAC code 
HMehani et aL][2007l [van der Hoist et al.|[2008 l) with 
Synge-type equation of state. We take a base resolu- 
tion of 120 X 120, allow for 5 grid levels, reaching an 
effective resolution of 1920^. Shorter timescale runs at 
even higher resolutions were done to confirm that the 
dominant initial large scale structure development is 
adequately resolved, although more fine scale features 
inevitably turn up. The simulations are performed typ- 
ically using 120 processors for about 2 days per case. 
We add some white noise both at the initial time and 
at time t - \ when various waves have already devel- 
oped (this latter addition may not be essential to the 
evolution). The fact that we use cartesian coordinates 
could give preference to instabilities with mode num- 
ber proportional lo m - A character (the case studied 
in Meliani & Keppens ( 2007 ) was therefore confirmed 



separately in cylindrical coordinates, with the overall 
mode number dominance of m - A recovered). How- 
ever, in all cases we study here, the instabilities dom- 
inating the dynamics and evolution of the jet have a 
clear mixture of many mode numbers. 

3. Results 

In the actual simulations performed in this work, 
the magnetic field in the jet is taken purely poloidal. 
As explained before, our five simulations differ in their 
magnetic configuration: (A) is hydrodynamic; (Bl) 
and (B2) has only the inner jet magnetized; (C) has 
only the outer jet magnetized; while in (D), both out- 
flows are magnetized. This represents the main differ- 
ence between the five cases, and while all cases will 
develop fairly complex nonlinear evolutions governed 
by multiple, interacting instabilities, our main aim is to 
determine which configuration can result in a clearly 
sustained two-component jet flow over a sufficiently 
long distance (time). It will turn out that cases (A), 
(C), and (D) all result in deceleration due to mixing 
between the two components, leading to decollimation 
of the jet, while only cases (Bl) and (B2) convincingly 
maintain their two-component character. For later ref- 
erence, important models parameters are listed in Ta- 
ble [T] We now continue to discuss the complex non- 
linear evolution of the five cases in some detail. 
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Table 1 : The most relevant characteristics and parameters for all models investigated. In addition to these 
tabulated values, the number density of the inner jet is fixed at 6.92 x lO^^cm while Vz = 0.99939c 
and = 0.01c. Constant values for the outer jet are number density 119.94cm"^, with - 0.9428c and 
= 10"^c. The external medium is always static and unmagnetized, and has number density lO^^cm"^. 







inner jet 


outer jet 


effective inertia ratio 


case 




Bz,m(0.1375Gauss) 




y^ph + B? ratio out/in 


A 


2.3 


0.0 


0.0 


0.15 


Bl 


0.1 


0.789 


0.0 


3.2 


B2 


10-4 


2.28 


0.0 


18.3 


C 


2.3 


0.0 


2.323 


0.15 


D 


2 


0.789 


2.323 


0.17 



3.1. Case (A) 

According to the Rayleigh criteria mentioned, the 
interface between the two components is always unsta- 
ble. Also, the effective inertia of the inner jet is higher 
than in the outer jet. In fact, at the interface R - R^, 
we have initially (y^ph^ x 0.07 (y^p h^. . In Ta- 
ble [T] a f = ratio is given using a value midway the 
outer jet divided by the axial value. The overall effec- 
tive inertia contrast is very different from the case we 



investigated in Mehani & Keppens ( 2007j i , where the 
ratio between (mean value of) effective inertia of the 
inner component jet to the effective inertia of the outer 
component jet was such that (y^pli^ ^ » {j^ph)- ■ 

Initially, a linear surface mode develops at this in- 
terface R - Ri„. As a direct consequence, a small, ra- 
dially extended shear region with low effective inertia 
y^p h forms at this location. At both bounding inter- 
faces of this small shear layer, one with the outer and 
one with the inner jet, small wavelength instabilities 
develop. Meanwhile, acoustic waves from the inner 
interface propagate inwards, reflect on the jet axis, and 
then give rise to a Kelvin-Helmholtz-like body mode 
with dominant azimuthal mode number m ~ 8 in the 
inner jet. Meanwhile, the surface instability at the 
interface between inner jet and shear shell reaches a 
clearly non-linear phase. The overall interaction then 
gives rise to the growth of a relativistically enhanced 
Rayleigh-Taylor-type instability, having a body mode 
character propagating into the inner jet (Fig. |2] top). 
This instability regime shows about four arms at first, 
that collapse to three arms which progressively propa- 
gate inwards. As will be detailed later on, this instabil- 
ity can be explained from the fact that in this case (A), 
fluid elements in the shear layer have a lower effec- 
tive inertia y^p h and higher pressure than pre-existing 
material in the inner jet. 



We could qualitatively describe this relativistic 
Rayleigh-Taylor instability as follows. A fluid ele- 
ment in the shear layer is at first typically rotating with 
a speed of order of ~ 0.01, the initial rotation speed 
at the inbound interface. The lower inertia of the shear 
shell relative to the inner jet makes that transfer of an- 
gular momentum from the inner wind/jet region to the 
shear layer is rather efficient, causing a fast rotational 
speedup of the shear layer matter. At the same time, 
the centrifugal force acting on a fluid element moving 
inwards varies with IjR. Because the fluid inertia in 
the shear region ends up at about a twelve times lower 
value than the inertia in the inner jet, the centrifugal 
force acting on this fluid element is much weaker than 
in the inner wind. At the same time, the pressure 
in the shear layer fluid element is higher than in the 
inner wind. Therefore, the element will expand in- 
wards until the pressure gradient gets balanced by the 
centrifugal force. In this particular simulation, three 
prominent fingers form and start propagating inwards. 
However, as the initial pressure profile decreases to- 
wards the jet axis, the centrifugal force in the inner jet 
decreases inward. With such a configuration, we find 
that in the phase where the intruding, rotating fingers 
of the Rayleigh-Taylor mode propagate inward, they 
will nearly reach the jet axis. Very close to the axis, 
the centrifugal force acting on them gets compensated 
by the pressure force. In the phase that follows, the ro- 
tating fingers get deflected sideways in the direction of 
rotation. After less than half a full rotation of the inner 
jet (at t ~ 32.6year, see Fig. |2] top), this relativistic 
Rayleigh-Taylor mode then dominates the body mode 
instability in the inner jet. The Rayleigh-Taylor fin- 
gers propagating inwards also compress the inner jet. 
Since angular momentum is conserved, the centrifugal 
force increases and the inner jet fluid moves outwards 
towards the shear shell. This means that we get en- 
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Fig. 2. — Case (A), the purely hydrodynamical jet, 
showing logarithm of proper density at times (top) t - 
32.6year: developement of relativistically enhanced 
Rayleigh-Taylor-type instability propagating inward in 
the inner jet, (center) t - 65.3year: merging of the 
three Rayleigh-Taylor fingers at the jet axis and (down) 
t - 163year: decollimation of the jet 

(one fuU rotation of the inner jet corresponds to 
t = 65.3year). 



hanced angular momentum transfer from the inner jet 
to the shear layer shell, which accelerates the rota- 
tional speed of the shear shell. We find that after less 
than one total rotation of the inner jet (at t ~ 65.3year, 
see Fig. |2] center), the three Rayleigh-Taylor fingers 
nearly merge, making the inner jet appearance totally 
dominated by the (growing) shear shell. 

An analysis of this relativistic Rayleigh-Taylor-type 
instability is described seperately in Section |4] The 
low effective inertia of fluid in the shear region makes 
the interface between the shear flow and the inner jet 
unstable to this instability, where centrifugal forces 
act as an effective (radially outwards pointing) gravity. 
Once initiated, these grow faster than the previously 
formed Kelvin-Helmholtz body mode. A number of 
'arms' develop from the shear region, having lower an- 
gular momentum than the fluid from the inner region, 
these propagate inwards, while inner jet fluid having 
higher angular momentum propagates outwards. The 
end result is transfer of energy and angular momen- 
tum from the inner jet to the growing shear region. As 
a consequence of this complex mixing of jet and shear 
layer matter, we find significant deceleration of the in- 
ner jet, while the density and the radius of the inner jet 
plus shear region increases (this is quantified later in 
Fig. [8). 

The external boundary of the growing shear region 
at R ^ Rin which borders the outer cold jet is also 
Kelvin-Helmholtz unstable. For the case at hand, in a 
first phase, a larger scale Kelvin-Helmholtz type mode 
develops with 8 undulations initially, eventually con- 
verging to form about 5 large scale structures (see 
Fig.|2] center). This is different from the case studied 
before in |Meliani & Keppens ( 2007 | l, where the most 
prominent feature in the non-linear evolution was the 
development of a Kelvin-Helmholtz type mode with 4 
arms. This mode transfers angular momentum from 
the shear region to the outer jet and decollimates the 
jet. At smaller spatial scale, but also at the interface 
between the shear region and the outer jet, other in- 
stabilities develop as well. These are again more cen- 
trifugally driven, as angular momentum and pressure 
increase in the shear region. Thus also at this interface, 
the centrifugal force becomes locally higher than in the 
outer outflow. Small bubble-like protrusions from the 
shear shell enter the outer outflow. However, since the 
pressure of the outer jet slowly increases outward, the 
growth of this smaller scale instability is stopped when 
the centrifugal force of the protruding bubbles gets bal- 
anced by the pressure of the outer jet. Therefore, this 



instability remains weak and the large-scale dynam- 
ics of the interaction between the shear layer and the 
outer jet remains dominated by the Kelvin-Helmholtz 
instability. The small scale instabilities do influence 
the dynamics, by increasing the efficiency of the an- 
gular momentum transfer from the Kelvin-Helmholtz 
structures to the outer jet. Finally, also the outer inter- 
face at R X Raut between outer jet and the hot, dilute 
medium is Rayleigh-Taylor unstable, and is dominated 
by small wavelength perturbations. 

As the overall outcome of all the interacting in- 
stabilities, the two-component jet decollimates, such 
that the jet radius increases to R^ut ~ 0.28pc (see 
Fig. |9]l. The inner component jet also spreads and its 
radius reaches R ~ O.lSpc (see Fig. [8]). Also dur- 
ing the simulated 163 years in physical time (Fig. |2] 
down), the rotation speed of both components de- 
creases and the interaction between the shear region 
and the outer component is dominated by Kelvin- 
Helmholtz induced behaviour, together with the rela- 
tivistically enhanced Rayleigh-Taylor instability of the 
inner rotating fluid. The main result of this interaction 
between the two components is the dramatic decelera- 
tion of the inner jet. In fact, the relativistic Rayleigh- 
Taylor-type instability leads to a deceleration of the 
inner component where the Lorentz factor of the in- 
ner component drops to y ~ 8 (quantified later for all 
models in Fig.|7]i. 

3.2. Cases (Bl) and (B2) 

We now describe the two-component jet evolution 
with a magnetized inner component. The effective 
inertia in the inner jet is lower than in the outer jet 
with initial contrasts (y^ph^ ^ x 3.2 (j^pli + B^). in 

case (B 1 ) and (j^p ~ 1 8 (y^p h + B^. in the case 
(B2), due to the initial lower pressure at the jet axis and 
the initial different distribution over thermal and mag- 
netic energy in them. These cases are then such that 
the inner interface is now stable against the dominat- 
ing Rayleigh-Taylor type instability described above 
and explained analytically in Section]?] Still, the in- 
terface at R ^ Rin between the inner and outer jet is 
subject to Kelvin-Helmholtz instability because of the 
differential rotation. This modifies the shape of the in- 
terface surface, leading to non-axisymmetric reflection 
of magnetohydrodynamic waves propagating through 
the inner jet. These waves in turn disturb locally the 
initial equilibrium between total pressure gradient and 
centrifugal force. 




Fig. 3. — Case (Bl), only the inner jet is magnetized. 
Logarithm of density at (top) t = 32.6year: devel- 
opement of Kelvin-Helmholtz body mode instability 
at the interface between the inner and outer jet, (cen- 
ter) t = 65.3year: inward growing of the shear shell in 
the inner jet and (down) t - 163year: developement 
of Rayleigh-Taylor at the outer jet interface, with one 
rotation of the inner jet completed att- 65.3year. 
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Time: 32.6 (year) 




Time: 65.2 (year) 




Time: 163. (year) 




Fig. 4. — Case (B2), only the inner jet is magnetized. 
Logarithm of density at (top) t = 32.6year: devel- 
opement of Kelvin-Helmholtz body mode instability 
at the interface between the inner and outer jet, (cen- 
ter) t = 65.3year and (down) t = 163year: the inner jet 
component is deplaced from its on-axis position due to 
significant non-axisymmetric mode development, with 
one rotation of the inner jet completed at / = 65.3year. 



During the evolution, some spikes develop at the in- 
terface, where the outer jet locally interchanges with 
the inner jet. These spikes are accelerated in the 
toroidal direction by the faster rotating inner jet. The 
centrifugal force acting on these slows down their in- 
ward expansion, and they then mainly propagate in the 
toroidal direction. Their interaction with the inner jet 
material induces Kelvin-Helmholtz body mode insta- 
bility, with a spiral pattern forming, this time having 3 
arms in the two cases (see Fig.|3] top and Fig.|4] top). 
They extract some angular momentum from the outer 
regions of the inner jet, slowing its rotation. Again, 
a shear shell then forms in this region, continuously 
extracting angular momentum from the inner jet/wind. 
The spikes then slowly expand to the jet axis. The in- 
ner wind gets compressed by the shear layer How- 
ever, these mainly Kelvin-Helmholtz type instabilities 
turn out to have lower efficiency in extracting angular 
momentum than the relativistic Rayleigh-Taylor insta- 
bility encountered in case (A). The compression of the 
inner jet by the shear region is followed by a modest 
expansion. In the case (B2), the extraction of the angu- 
lar momentum is weaker than in the case (Bl), because 
in case (B2) the inner jet has lower inertia. Mean- 
while, the angular momentum extracted from the inner 
jet leads to an outwards extension of the shear shell 
that formed at Ri^. At the interface between this grow- 
ing shear shell and the outer jet, a Kelvin-Helmholtz 
instability surface mode with azimuthal mode num- 
ber m ~ 3 character then develops. Also at this in- 
terface, smaller scale Rayleigh-Taylor instability pro- 
ceeds. The shear shell has an effective inertia which is 
lower than the outer jet inertia. As locally centrifugal 
forces dominate, low inertia 'bubbles' radially extend 
toward the outer jet. The growth rate of these bubbles 
is rather slow, since the pressure of the outer jet in- 
creases mildly outwards and the slowly rotating outer 
jet breaks their rotation. 

At the outer surface /?out of the outer jet in the case 
(Bl) (see Fig.[3]l, we also witness Rayleigh-Taylor in- 
stability, as in case (A). On the other hand, case (B2) 
shows little evidence for such instability at the outer jet 
surface (see Fig.|4]l, which is due to the different exter- 
nal pressure conditions there. Throughout the simula- 
tion, the jet in (Bl) and (B2) remains collimated by the 
outer component, which compresses the inner compo- 
nent and shear region. The inner component jet radius 
in case (Bl) at 163 years remains lower than 0.07pc 
and in case (B2) even below about 0.04pc (see Fig.[8]l. 
However, in case (Bl) the Rayleigh-Taylor instability 
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at the surface of the outer region forms an extended 
sheat, which gives an increase to the apparant total jet 
radius 7?jet = 0.2pc (see Fig.|9jl, whereas the jet in case 
(B2) does not show any sign of decollimation. 

Both the inner jet and the shear region end up mag- 
netized. The main difference between these cases and 
all others studied here, is that despite all small-scale 
instability development, the inner jet decelerates lit- 
tle, dropping to an average Lorentz factor of about 20 
(see Fig. |7]l. The two-component jets for case (Bl) 
and (B2) remain clearly separable in inner and outer jet 
components (see Figs. [3] down-|4j down), which is not 
the case for all other evolutions shown here. In both 
cases (Bl) and (B2), we do find that the inner jet com- 
ponent is deplaced from its on-axis position due do 
significant non-axisymmetric mode development. The 
jet stratification converges to a structure with an inner 
fast, magnetized spine having a Lorentz factor of about 
20. The spine is surrounded by a shear shell being 100 
denser than the spine and with lower Lorentz factor 
and low magnetisation. The difference between the 
two cases is that the high magnetic pressure and low 
thermal pressure in the inner jet for case (B2) increases 
its non-linear stability as compared to case (Bl). 

3.3. Case (C) 

In this case, only the outer jet is magnetized with a 
purely poloidal magnetic field. As mentioned before, 
the rotation of the jet is subsonic, such that acoustic 
waves which start to propagate from the inner-outer 
jet interface at different stages in the evolution propa- 
gate inwards and are reflected about 200 times during 
one rotation period of the inner jet. Our simulations 
follow about 2.5 full inner jet rotations. As the inter- 
face between the two components is unstable to vari- 
ous mode types as already encountered in the previous 
cases, its shape changes and the acoustic wave prop- 
agation becomes strongly non-axisymmetric. These 
waves play a clear role in disturbing the initial axisym- 
metry of the inner jet, and mitigate the development of 
non-axisymmetric body mode instabilities in the inner 
jet. 

In this case (C), in a first stage, linear small- 
wavelength (compared to case (A)) Kelvin-Helmholtz 
instabilities develop at the interface R R^„ between 
the two components. These reach a non-linear phase 
in less than a tenth of a full rotation. Spikes of the 
outer jet, which are this time magnetized, arise in the 
inner jet. The large ratio of the fluid inertia between 



Time: 32.6 (year) 




Fig. 5. — Case (C), only the outer jet is magnetized. 
Logarithm of density at (top) t - 32.6year: Rayleigh- 
Taylor instability develops at interface between the 
two component inducing rotating 'bubbles' that pro- 
trude inward, (center) t = 65.3year: 'bubbles' merge 
at the jet axis and inner jet decelerates and becomes 
turbulent and (down) t = 163year: The shear shell ex- 
pand outward in the outer jet and decoUimates the jet, 
with one full inner jet rotation at f = 65.3year. 
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the outer and inner jet makes these spikes once more 
unstable to the centrifugally mediated, relativistically 
enhanced Rayleigh-Taylor instability seen previously 
in case (A). This gives rise to the formation of rotating 
'bubbles' that protrude inward. In the case at hand, 
about ten such bubbles form (see Fig. [5] top), hav- 
ing smaller size compared to those seen to emerge in 
case (A). This difference results due to the difference 
in the Kelvin-Helmholtz surface mode development 
mentioned earlier The initial smaller size of the pro- 
truding bubbles makes their growth rate slower than in 
the purely hydrodynamical case (A), until they reach 
the merging phase where in this case about three larger 
size bubbles form and converge to the jet axis. Mean- 
while, these bubbles compress and push fluid elements 
from the inner jet to a shear layer (see Fig. |5j;enter). 
Thus once again, the larger-scale mode structures act 
to transfer angular momentum from the inner jet to the 
shear region. 

Since the bubbles which move radially inwards 
form at the shear shell, their ratio of magnetic pressure 
to thermal pressure is lower 

than in the outer jet. In fact, the shear region is made 
up of fluid from the inner and the outer jet, and its ther- 
mal pressure dominates total pressure. The external 
surface of the shear shell is Kelvin-Helmholtz unstable 
and about 6 arms form (see Fig.|5] center). From each 
one, a centrifugally driven spike propagates outward. 
In fact, angular momentum extracted from the inner 
jet in the end gets transferred to the external part of the 
shear region. Then, locally the centrifugal force at the 
interface between the shear shell and the outer jet acts 
to develop spikes which decollimate the jet (see Fig.|5] 
down). The total jet radius increases from O.lpc at the 
initial stage to 0.189pc (Fig.|9|, after the jet has prop- 
agated for a distance of 50pc. In the same time, the 
inner jet spreads as well, decelerates and its Lorentz 
factor drops to about 8. At the endstate, the jet is con- 
stituted of an inner turbulent component with a radius 
of O.lSpc (see Fig. [8]) and a small outer component. 
The final rotation speed profile increases from the axis 
to 0.006 at 0.04pc, further decreasing outwards. 

3.4. Case (D) 

In this last case, both inner and outer jet are magne- 
tized. Initially, from the discontinuity between the two 
components, fast magnetosonic waves propagate both 
outwards into the outer jet, and inwards into the inner 
jet. The end result of the initial readjustment is that a 
small shear region with lower magnetization than both 




Fig. 6. — Case (D), both jet components magnetized. 
Logarithm of density at (top) f - 32.6year: 13 bubbles 
arise, moving inwards into the inner jet, (center) t = 
65.3year: Rayleigh-Taylor bubbles emerge at the jet 
axis and the inner jet becomes turbulent and (down) 
t = 163year: The shear region grows and changes the 
full jet shape to elliptic, with one full rotation of the 
inner jet corresponding to t - 65.3year. 
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The evolution of mean Lorentz factor 



The evokition of the inner jet component radius 




t (year) 



Fig. 7. — Time evolution of the mean Lorentz factor of 
the inner jet in the five cases. For the cases (A), (C), 
and (D) we distinguish between the inner and outer jet 
according to the polytropic index and Lorentz factor 
In the cases (Bl) and (B2), we use the same condition 
for the Lorentz factor to distinguish between the inner 
and outer jet, as well as the extra condition that the 
inner jet is magnetized. 



component develops at the interface R « Also 
this configuration is Kelvin-Helmholtz unstable. At 
the surface from the shear region roughly 13 bubbles 
arise, moving inwards into the inner jet (see Fig. [6] 
top). They merge to form 3 large arms expanding to- 
wards the jet axis (see Fig. |6] center). Fluid from the 
inner jet is pushed outwards, thereby extracting angu- 
lar momentum from the inner jet. Since the structures 
form at the shear region, the magnetic pressure inside 
the bubbles is lower than in the inner and outer com- 
ponent. However, the thermal pressure is increased 
in the shear region. During the evolution, dominant 
forces are magnetic pressure and thermal pressure that 
balance each other, while the contribution of the cen- 
trifugal force is low. However, the centrifugal force 
remains responsible for triggering instabilities. As in 
case (C), the shear region grows with eventual almost 
complete destruction of the inner jet component (see 
Fig. |6] center). In the last phase simulated, the dom- 
inating inner shear region becomes totally turbulent, 
and strong vortices form, changing the full jet shape 
to elliptic (Fig. |6j down). At this phase the average 
Lorentz factor of the inner jet has dropped down to 
Tmax ~ 7 (see Fig.|7]i. 



Fig. 8. — Time evolution of the effective inner jet ra- 
dius in all cases. For the cases (A), (C), and (D) we 
distinguish between the inner and outer jet according 
to the polytropic index and Lorentz factor. For the 
cases (Bl) and (B2) we use the same condition for 
the Lorentz factor and the condition that the inner and 
shear component are magnetized. The eff'ective inner 
jet radius is yfS^Jn where 5 in is the surface thus found 
for the inner jet. The jets (Bl, B2) with low inner ki- 
netic energy flux contribution are more stable and re- 
main relativistic and the jets (A, C, D) with high inner 
kinetic energy flux contribution are unstable and de- 
celerate 



4. Relativistically enhanced Rayleigh-Taylor 

In cases (A), (C), and (D), a relativistic, Rayleigh- 
Taylor type instability developed. Here, we first de- 
scribe in a qualitative way how this instability ensues, 
and then continue to explain it using an approximate 
linearization. 

At the interface between the two outflow compo- 
nents, we mentioned that (linear and then non-linear) 
surface modes develop, due to velocity shear hence 
of Kelvin-Helmholtz type. Fluid of lower inertia (low 
y^ph) and higher pressure then penetrates the inner jet. 
Since the inertia in the penetrating spikes is lower than 
in the surrounding inner jet material, the centrifugal 
force acting on these is weaker. The spikes then ex- 
pand inwards until the pressure is balanced by the cen- 
trifugal force. They expand sideways as well, again on 
the basis of pressure balance and ram pressure argu- 
ments. Due to the big difference in inertia, transfer of 
angular momentum from the inner jet to these inward 
penetrating spikes is very efficient. Rather quickly 
then, theii- rotation speed reaches the rotation speed of 
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the inner jet at the initial interface, set by v^ in- Dur- 
ing its further inward expansion, this rotation speed re- 
mains constant. However, the low inertia of the spikes 
keeps their angular momentum lower than the inner jet 
matter. 

To explain this phenomenon semi-analytically, we 
confine our attention to a perturbation depending only 
on the radial direction R. Let dR = v^dt represent 
the radial displacement of the contact surface between 
the two components during dt, and vr the perturbation 
speed. The momentum equation governing the fluid 
near our initial equilibrium writes as 



variation of the total pressure can be argued to lead to 



d_ 
dt 



+ y-v y+V;,.o.3i+V^+... = 0, 
ot 

(14) 

with ptotM - p+\- In this equation we already ignored 
the lab frame contribution of the charge separation, 
which is justified for rotational flows far within the 
hght cylinder, as is the case in our model here [^More- 
over, the magnetic topology taken here does not in- 
volve magnetic pinching, hence magnetic effects only 
add to total pressure and to the effective inertia, where 
also « Vz is used. Then the radial component 
of Eq. ( 14 1 can be further approximated, if we ne- 



glect temporal variation of the total pressure and use 



dv] 



^=0, to get 



dt 



R 



dp, 



total 



dR 



(15) 



According to the initial condition, the inertia ph 
of the fluid varies slowly in the inner and outer jet. 
We can thus argue that both inertia and toroidal speed 
in each fluid element varies slowly when it undergoes 
a small radial displacement ^ - Ri^ - R - J v^dt = 
- J VRdt from the initial interface position Ri^, with 
I ^ |« Ri„. Using the main equilibrium balance be- 
tween centrifugal force and total pressure gradient, the 



'Note that we mentioned a jet radius of O.lpc at a distance of Ipc for 
M87. Tlie estimated light cylinder radius for M87 is smaller than 
O.lpc if a large fraction of the angular momentum (then Poynting 
flux) extracted from the accretion by the magnetic field is transferred 
to the stream lines during the acceleration phase by the ideal RMHD 
mechanisms ^Begelman & Li|1994| . But in our model we assume 
that the jet is accelerated also by non-ideal RMHD mechanisms with 
decay of the toroidal magnetic field by dissipation/reconnection, 
leading to a fast central spine jet, and overall low rotation profiles 
(weak transfer of angular momentum from magnetic fields to stream 
lines). 



d J Vfd^ 



dt 



Rm 



(16) 



To get an approximate dispersion relation, we assume 
that to first order, the perturbation speed is potential, 
i.e. = (VT)^ and that we have ^' cc exp(At-k \ ( |). 
Noting that the displacement C - J ^'('^^ must be iden- 
tical for inner/outer regions, and using total pressure 
and displacement continuity arguments, we then get 
the essential proportionality relation between the in- 
stability growth rate and wave number 

A' cc k [{r'ph + Bl\^ - [y^ph + B^;)J . (17) 

This approximative dispersion relation indicates that 
for two-component stability, we need < Q rt- 

quiring diat (y^p/z + B^) > [j^ph + B^.. This 
confirms that the interface between the two rotating 
components is stable against a centrifugally driven, 
relativistic Rayleigh-Taylor instability, when the ef- 
fective inertia of the outer component is higher than 
the effective inertia of the inner component. This 
important result correlates well with the results of 
all simulations in section [3] as well as in our pre- 
vious publication (Meliani & Keppens 2007 1. In 
fact, this equation explains why cases (A), (C), (D) 
are relativistically Rayleigh-Taylor unstable since we 
find typical contrasts during the evolution in (A) 
(yV/^ + fiz)„„, ~ 0.75(r-p/^ + fiz)i„' in (Q and (D) 
{y^ ph + B^ K Q.l{y^ ph + B^. , whereas cases 
(Bl) and (B2) are relativistically Rayleigh-Taylor sta- 
ble since [y^ ph + B^)^^^ ~ lQ{y^ph + B^).^^ in (Bl) 

and {y^ph + Bl) ~ 20{y^ ph + B^). in (B2). Initial 
eff'ective inertia ratios derived from the initial condi- 
tions are mentioned in Table [T] and it is clear that we 
can set up stable versus unstable cases by varying the 
contrast of effective inertia between inner versus outer 
jet. This also works in pure hydro, and indeed the pure 
hydro case from [Meliani & Keppens| ( f2007| does not 
suffer from this newly discovered instability. For the 
purely poloidal field configurations studied here, mag- 
netic field effects are prominent in total pressure and 
effective inertia alone. Note finally that it is a truly 
relativistic effect, since the same argument in classi- 
cal MHD just involves the density difference between 
outer versus inner jet. Under the conditions of a light 
inner jet with heavy outer jet taken here, a classical 
variant of relation (17i involves only pout - Pin and 



would predict stability. 
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The evolution of the jet radius 



Case A 




Case B1 + 




Case B2 x 




CaseC 




CaseD 













t (year) 

Fig. 9. — Time evolution of the effective jet radius in 
all cases. The effective jet radius is ^JS]^J7: where 5 jet 
is the surface occupied by the full (inner+outer) jet. 
The jets (A, C, D) with high inner kinetic energy flux 
contribution decollimate. 



5. Discussion and conclusions 

We examined five configurations of magnetized 
two-component jets. All share the same density ra- 
tio between inner and outer component and identical 
rotation profiles. The magnetic and thermal pressure 
configuration in each model differs, though. This in 
turn translates to different distributions of the (total, 
fixed) kinetic energy flux over the inner and outer jet 
component. In fact, in case (A) the kinetic energy 
flux in the inner component is about 10%, while it is 
0.7% for case (Bl) and 0.5% in case (B2), and around 
38% for cases (C) and (D). For the outer component, 
we thus have 90% in case (A), 99.3% in case (Bl) 
and 99.5% in case (B2), and only 62% for cases (C) 
and (D). The most important difference between the 
two model categories is then: cases (A, C, D) have 
an inner jet component with higher inertia ph + 
than their outer jet component, while cases (Bl) and 
(B2) have an inner jet component inertia which ends 
up lower than in the outer jet component. From the 
detailed analysis of the simulations, as well as from 
the approximate stability analysis, this criterion distin- 
guishes between cases where relativistically enhanced 
Rayleigh-Taylor modes ultimately dominate the evolu- 
tion, leading to complete mixing of both components 
and inner jet deceleration. This is quantified most 
clearly by showing the time evolution of the mean 



Lorentz factor over the inner jet region for all cases 
in Fig.|7] This requires a clear criterion to distinguish 
inner versus outer jets in the turbulent evolutions. In 
cases (A), (C) and (D), we locate the outer jet com- 
ponent as having a Lorentz factor 2.5 < y < 3.5 and 
effective polytropic index Fefr > 3/2. The inner com- 
ponent is the region defined by a Lorentz factor y > 3.5 
and effective polytropic index Feff < 3/2. In the cases 
(Bl) and (B2), the effective polytropic index in the 
inner and outer jet can be locally of the same order 
there, hence we use that the inner jet is magnetized. In 
the cases (Bl) and (B2), the inner component jet and 
shear region are not totally mixed during the evolu- 
tion, since the inner component is compressed and has 
higher magnetization and Lorentz factor. Thus to dis- 
tinguish the inner component jet we put a condition on 
magnetic field strength > Bz,initiai/2, where Bz.imtiai 
is the magnetic strength assumed initially. Under these 
precise quantifications of inner/outer jet regions. Fig. [7] 
demonstrates clearly how stable cases (for the rela- 
tivistically enhanced Rayleigh-Taylor modes) remain 
at high speed, while unstable cases decelerate. Using 
the same means to distinguish inner versus outer jet 
regions at all times, we can quantify the inner jet ra- 
dius for all cases, as well as the total jet radius for all 
cases. These are shown in Figs. 8]|9 and quantify the 
decollimation effects discussed in Section [3] 

During the entire evolutions, the toroidal and radial 
speeds remain weak as we have typically a maximal 
Vr < 0.01 and < 0.04. This means that the conti'i- 
bution of the laboratory frame charge separation force 
PeE to the Lorentz force is negligible at all times. Un- 
der these conditions, in both radial and toroidal direc- 
tions, the contribution of magnetic energy to fluid in- 
ertia is 1/y^ weaker than the contribution of the mag- 
netic pressure to the total pressure. This explains why 
cases (Bl) and (B2) are then more stable than the other 
cases. Despite the fact that case (B2) has similar ax- 
ial total pressure than other cases (A, C, D), the low 
contribution of thermal energy to total pressure makes 
the effective inertia ratio between the inner and outer 
jet low. This two-component jet is then stable against 
the relativistically enhanced Rayleigh-Taylor instabil- 
ity. Extraction of angular momentum and energy from 
the inner jet to the shear shell is less efficient in cases 
(Bl) and (B2) than in all other cases. As a result, in 
cases (Bl) and (B2) the inner relativistic jet with high 
Lorentz factor yin ~ 20 persists. In the other cases, the 
kinetic energy flux in the inner jet was initially rela- 
tively high, making them unstable and leading to de- 
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celeration to Lorentz factors around 8. 

We investigated the stability of two-component jets 
beyond the launching region, where both components 
are rotating differently and a clear two-component 
structure exists. We initialized this model in accord 
with magneto-centrifugal models for jet generation, 
also using the observed analogy between radio source 
jets and two-component jets in young stellar objects, 
where the rotation within the jet can actually be ob- 
served. We performed five very high resolution simu- 
lations of magnetized two-component jets with various 
magnetization and kinetic energy flux stratifications. 
The two-component jets with a low inner kinetic en- 
ergy flux contribution are more stable and remain rel- 
ativistic for long distances, whereas jets with a highly 
contributing inner jet to the total jet kinetic energy flux, 
are subject to a relativistic Rayleigh-Taylor type insta- 
bility. This instability turns out to be very efficient to 
decollimate and decelerate the inner jet. Jets that are 
subject to this instability become turbulent after prop- 
agating for a distance of about 30pc. 

This new result on two-component jet models is 
important because it can explain the classification of 
radio sources in Fanaroff-Riley 1/11 categories accord- 
ing the energy stratification of the inner jet. This ulti- 
mately relates to the jet launch region and the prop- 
erties of the inner accretion disk. In fact, by anal- 
ogy between the FRI/FRII classification and the results 
of our model, an FRI jet would correspond to a two- 
component jet with a high energy flux contribution 
from the inner jet, whereas the FRII jet corresponds to 
relatively low energy fluxes in the inner jet. The model 
we propose here to explain the FRl/FRJl dichotomy is 



diff'erent from the model we proposed earlier ( Meliani 



|et al.|20d8 ^ where the transition occurs due to external 
density stratification. That model explains the group 
of peculiar "HYbrid MOrphology Radio Source" (HY- 



MORS) ( |Gopal-Krishna & W iita 2000) which appear 
to have an PR II type on one side and an FR 1 type dif- 
fuse radio lobe on the other side of the active nucleus. 
Since the launch conditions on each side are presum- 
ably similar in these kind of radio sources, the differ- 
ent Fanaroff-Riley morphologies on either side must 
be attributed to the change in the properties of the am- 



bient media, as shown convincingly in ,Meliani et al. 
( 2008[ l. The results of the present paper nicely com- 



plement these earlier findings with a quantifiable role 
of the central engine contribution. 

We currently continue this study of the interaction 
between two component jets in full 3D. We thereby 



intend to explore the relative influence of azimuthal 
versus longitudinal instabilities for realistic multi- 
component jets. Another extension is to allow for 
aximuthal magnetic fields, in accord with the initial 
profiles as given generally in this paper. It then rre- 
mains to be shown that (1) the newly discovered insta- 
bility persists in 3D hydro and magnetohydrodynamic 
configurations, where the potential role of axial mode 
development is incorporated, and introduces helical 
jet axis displacements; (2) how the instability gets 
modified (stabilized or destabilized) by the inclusion 
of toroidal field components, first in 2.5D neglecting 
helical axis displacements, and consecutively in 3D, 
where also current-driven kinks may occur. 
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